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Abstract 

Unexpected accelerator modes were recently observed experimentally for cold cesium 
atoms when driven in the presence of gravity. A detailed theoretical explanation of this 
quantum effect is presented here. The theory makes use of invariance properties of the 
system, that are similar to the ones of solids, leading to a separation into independent 
kicked rotor problems. The analytical solution makes use of a limit that is very similar 
to the semiclassical limit, but with a small parameter that is not Planck's constant, but 
rather the detuning from the frequency that is resonant in absence of gravity. 



1 Introduction 

The kicked rotor is a standard system used in the investigation of classical Hamiltonian chaos 
and its manifestations in quantum mechanical systems @, |3[]. The motion is classically de- 
scribed by the Standard Map. The size of the chaotic component in the phase space of this map 
increases with the driving strength, and when the latter is sufficiently strong unbounded diffu- 
sion in action space takes place. Important deviations are however found for some values of the 
driving strength f|, [| f|, due to the onset of so-called accelerator modes, that produce linear, 
rather than diffusive, growth of momentum along orbits in a set of positive measure. Quanti- 
zation imposes remarkable modifications. For typical parameter values, the classical diffusion 
is quantally suppressed by a mechanism that is similar to Anderson localization in disordered 
solids U], ||. The accelerator modes decay in time due to quantum tunneling ||. Quantum res- 
onances are found when the natural frequency of the rotor is commensurate with the frequency 
of the driving |7]]. The quasienergy states are then extended in angular momentum, leading to 
ballistic (i.e., linear) growth of the latter in time. 
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The quantal suppression of classical diffusive transport first observed in the quantum kicked 
rotor is actually a more general phenomenon, now known as dynamical localization. Theoretical 
predictions |]] prompted the first experimental observations of this phenomenon for microwave 
driven Hydrogen and Hydrogen like atoms 0. 

The most direct experimental realization of the quantum kicked rotor is achieved in the field 



of atom optics, by a technique pioneered by Raizen and coworkers JII|. Laser-cooled Atoms 
(first Sodium and later Cesium) are driven by application of a standing electromagnetic wave. 
The frequency of the wave is slightly detuned from resonance, so a dipole moment is induced in 
the atom. This moment couples with the driving field, giving rise to a net force on the center 
of mass of the atoms, proportional to the square of the electric field [O, [T2J. As the wave is 
periodic in space, the atom is thus subjected to a periodic potential. The wave is turned on 
and off periodically in time, and the time it is on is much shorter then the time it is off. A 
realization of a periodically kicked particle is then obtained. In real experiments the duration 
of the kicks is always finite, setting a bound on the the momentum range wherein the 5-kicked 
model is applicable. At very large momentum the driving becomes adiabatic, leading to trivial 



classical and quantum localization in momentum [T_3|, [10| . 

The basic difference between the kicked particle and the kicked rotor is that the momentum 
of the particle is not discrete as is the angular momentum of the rotor. As we shall review in 
Sect. (2.2) , this difference is circumvented by Bloch theory. The spatial periodicity of the driving 
only allows for transitions between momenta that differ by integer multiples of hG, with 2ttG~ 1 
the spatial period of the driving potential. This implies conservation of quasi-momentum. The 
particle wavepacket is a continuous superposition of states with given quasi-momenta. The 
dynamics at any one fixed quasi-momentum is that of a rotor, which however differs from 
the standard kicked rotor because of a constant shift in the angular momentum eigenvalues, 
proportional to the given quasi-momentum. This modification of the kicked rotor dynamics is 
formally the same as that produced by a Aharonov-Bohm flux threading the rotor. It does not 



crucially affect dynamical localization ||14j| , so the latter carries over to the particle dynamics. In 
experiments JTU] it is found that for typical values of the parameters, an initial narrow gaussian 
distribution around zero momentum spreads into an exponential distribution, characteristic of 
dynamical or Anderson localization. 

The difference between the rotor's and the particle's dynamics due to the presence of a 
continuum of quasi-momenta is indeed crucial in what concerns quantum resonances at "com- 
mensurate" values of the kicking period, because these only occur at special values of quasi- 
momentum. Therefore, nearly all quasi-momenta involved in a particle's wave packet would not 
be in resonance [ I3| . In this paper we present an exact calculation, showing that the quadratic 
spread in momentum grows in this case linearly in time, in contrast to the situation found for 
the kicked rotor, where this spread is quadratic in time. A more detailed analysis of the kicked 



particle dynamics at resonance will be presented elsewhere [16 



In the above discussed experiments, gravity had but negligible effects, as driving of the 



atoms took place in the horizontal direction. In recent experiments |17|, 18, [19|| , that provide 
the subject of the theoretical analysis of the present paper, atoms were driven in the vertical 
direction, and gravity was found to produce remarkable effects. In the vicinity of the resonant 
frequencies of the kicked rotor, a new type of ballistic spread in momentum was experimentally 
observed [TR 05 1. A fraction of the atoms are steadily accelerated, at a rate which is 
faster or slower than the gravitational acceleration depending on what side of the resonance 



2 



the driving frequency is. Such atoms are exempt from the diffusive spread that takes place for 
the other atoms, and their acceleration depends on the difference between the driving and the 
resonant frequencies. [] In ref. [|T8| a physical explanation was given, and it was stressed that 
the phenomenon is resemblant of the accelerator modes in the Standard Map. The accelerating 
parts of the distributions were hence termed "quantum accelerator modes" , at once emphasizing 
that resemblance, and their purely quantal nature; in fact, they have no classical counterpart 
in the classical dynamics of the kicked particles in gravity.0 

In this paper we present a theory that explains the experimental results in terms of the 
exact quantum equations of motion, and allows for further predictions. The main result is 
that the quantum accelerator modes in presence of gravity do indeed correspond to accelerator 
modes of a classical map. This map is not, however, the one given by the proper classical limit 
(h — > 0). It emerges of a quasi-classical asymptotics, where the small parameter is not h but 
rather the detuning of the kicking frequency from the resonance of the kicked rotor. Though 
a formally simple variant of the Standard Map, it is endowed with a rich supply of accelerator 
modes and complex bifurcation patterns. 

Our analysis starts from the time-dependent Schrodinger equation for the kicked particle in 
the laboratory frame. A simple gauge transform reduces the corresponding Floquet operator to 
that of a particle kicked by a potential, that is quasi-periodic (and not just periodic) in space 
(Sect. (2.1)). The corresponding incommensuration parameter adds to the kicking period, in 
building a formidable mathematical problem. In particular, quasi-momentum is not conserved, 
preventing implementation of the Bloch theory. Moving to a free-falling frame removes this 
difficulty and restores decomposition into independent rotor problems, at the price of time- 
dependent Floquet operators (Sect. (2.3)). On such operators we work out the mentioned quasi- 
classical approximation at small detuning from resonance. 

Variants of the kicked-rotor, in which some parameter was allowed to change with time, 
have been considered earlier ||20|| , the issue being what degree of uncorrelatedness in the time- 
dependence is sufficient to destroy localization. For the quasi-periodic dependence of the present 
model, this issue is as yet unsolved, and depends on incommensuration properties. In a mathe- 
matical aside of this work, we prove that in the presence of gravity and at resonant values of the 
kicking period the particle energy (in the falling frame) grows linearly in time, in sufficiently 
incommensurate cases at least. 



2 Discrete-Time Quantum Dynamics 

The dynamics of the atoms that are falling as a result of gravity and are kicked by the external 
field is modeled by the time-dependent Hamiltonian: 

H(t) = 2^- MgX + kcos(GX) £ 6{t - tT), (1) 

t=— oo 

where t is the continuous time variable, the integer variable t counts the kicks, P, X are the 
momentum and the position operator respectively, M is the mass of the atom, 2ttG~ 1 is the 

1 The main experimental results are clearly presented in Figs 4 and 13 of Ref. [Q. 

2 For some parameter values, the latter dynamics does exhibit accelerator modes, which are however unrelated 
to the experimentally observed ones. 
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spatial period of the kicks, k is the kick strength, and T is the kicking period in time. The 
positive x— direction is that of the gravitational acceleration. Without changing notations, we 
rescale momentum P in units of hG, position X in units of G~ l , and mass in units of M. Then 
the energy E comes in units of h 2 G 2 /M, and time t in units of Mj {hG 2 ). The reduced Planck 
constant is equal to 1, and the Hamiltonian takes the following form: 

H(t) = — -X + k cos(X) V 6(t - tr) , (2) 

9 T * J 



where: 



k - - r - n - *W (3) 

In the above defined units, rj/r is the gravitational acceleration. The dynamics is fully charac- 
terized by the dimensionless parameters k,T,rj. 

In the following Dirac notations will be used: e.g., ip(x) = (x\ip) and ip(jp) = (p\ip) will 
denote the wave function in the position and in the momentum representation respectively. 

2.1 Floquet operators. 

The quantum evolution over a sequence of discrete times spaced by one period of the external 
periodic driving is obtained by repeated application of the Floquet operator. This is the single 
unitary operator which gives the evolution from any instant to the next one in the given sequence 
of times. Our discrete times will be the kicking times themselves Throughout the following, 
time is a discrete variable, given by the kick counter t. As the state discontinuously changes 
at the kicking times, we further specify the state at time t to be the one immediately after 
the t—th kick. The Floquet operator U is then found by integrating the Schrodinger equation 
(with the Hamiltonian (|2])) from t= + to t= r + : 

U = KF = e -*«*(*)p, (4) 

where K describes the kick, and F describes free fall inbetween kicks. With the present units, 



the energy eigenfunctions ue(p) = (p\E) of the particle in the gravity field read [fJTJ : 



/ T \ V2 3 

u e(p) = I - — J e ^ (jBp "^ r) . 
\27TriJ 

Therefore, apart from a constant phase factor, 

(p'\F\p") = JdE e- iEr u E (p')u* E (p") 
= 5{p'-p"-ri) e^W-W 

3 Other choices lead to different Floquet operators, which are nonetheless unitarily equivalent to the present 
one. 
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Moreover, 

oo 

(p\e- ikcos ^\p') = J2 K n 8{p-p'-n) 

n=—oo 

where K n = (—i) n J n {k) and J n (k) are the Bessel functions of the 1-st kind. Replacing in (||) 
we get the explicit form of the propagator in the laboratory frame: 

oo 

(Uif>)(p) = R n e-^^-a) 3 V(p - n - 77) . (5) 

n=— 00 

Note that K n = K_ n . A more transparent formulation is gained by introducing the operators: 

R = e- tTp2 ' 2 , S = e^' 2 

One may then write: 

U = SKRS = S^U'S , U' = S 2 KR (6) 

Thus U only differs by a unitary transformation (in fact a gauge transformation) from U', which 
has the simple form: 

(jl _ _ e i(vX-k cos(X)) e -irP 2 /2 ^ 

Thus formulated, the problem is that of a particle freely moving on a line, except for time- 
periodic kicks. The spatial dependence of the kicks is periodic when rj is rational, quasi-periodic 
otherwise. 



2.2 Quasi-momentum and Kicked Rotors. 

If 77 = 0, then ([7D is formally similar to the Kicked Rotor model [|lj], from which it however 
differs in one crucial respect: whereas the latter has the kicked particle moving on a circle, (0) 
has the particle moving on a line instead. A link between the two models is established by the 
spatial periodicity of the kicking potential. We review this well-known construction, because it 
plays a fundamental role in this paper. 

At rj = the evolution operator U commutes with spatial translations by multiples of 2tt. As 
is well known from Bloch theory, this enforces conservation of quasi-momentum. In our units, 
this is given by the fractional part of momentum, and will be denoted by (3. We then introduce 
a family of fictitious rotors (particles moving on a circle) parametrized by (3 £ [0, 1), with angle 
coordinate 9 (henceforth named (3— rotors), and denote their states. For integer n we 
denote \n) the angular momentum eigenstates of these rotors so that in the 9— representation 
(9\n) = (27r)~ 1//2 exp(m#). To states of the particle we associate states of the (3— rotors 
as follows: 

(9\*p) = _L^<n + /3h/>>e^. (8) 

V Z7T — 

n 

In the angular momentum representation, 

(nl^) = (n + m (9) 
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Note that |\Jjg) is not necessarily normalized to 1, even if 
particle is retrieved from the (3— rotor states via 



is. Conversely, the state of the 



(?#) 



^7T 



2jt 



d9 (6\Vp) e 



-ind 



n 



[p] ,p = {p} 



(10) 



where [p] , {p} denote the integer and the fractional part of the momentum p respectively. Using 
(|10| ) and Poisson's summation formula, one obtains the wave function in the x— representation: 



(x\ 



dpe ipx (p\ip) 



dpe^ x (x mod(27r))|* /3 ) 



Fixing a sharp value of (3 yields a spatially extended state (a Bloch wave ) for the particle, even 
with a normalizable rotor wave function. 



It follows from @ and @ (with r) — 0) that as evolves into U t \ 
into Up\fyp), where: 



U,3 = (KKp) , (n\JC\m) = K n 



(n\ilp\m) 



Eqn. fll2|) may also be written as follows: 



e -ifccos(e) e -ji(AT+/3) 2 



\E^) in turn evolves 
(12) 

(13) 



where Jv = Y2 n n \ n ) ( n l ^ s angular momentum operator: Jv = —%4a in the representation, 



At (3 = 0, eqn. ( p~3|) is the standard Floquet operator of the Kicked Rotor. At (3 ^ one obtains 
a variant of the Kicked Rotor, which has also been studied, (3 being typically given the physical 



meaning of an external magnetic flux||T4 



We have thus shown that in the presence of conserved quasi-momentum the particle dynam- 
ics can be determined as follows: given an initial particle (pure) state | -0 (0) ) one first computes 
the corresponding rotor states (^(O)) as above described. These separately evolve into \ty p{t)) 
at time t. The particle state \ip(t)) is finally reconstructed using (|10|) . The process of decompos- 
ing the particle dynamics in a bundle of rotors will henceforth be named "the Bloch- Wannier 
fibration" . It is the quantal counterpart of the classical process of "folding back" the particle 
trajectory onto a circle, by taking the x— coordinate mod(27r). 

The kinetic energy of the particle at time t is: 



£{t) 



- I dpp z \^{ P) t)[ 




1 oo 



dp (n + /3) 2 |^(n + /?,t)p 



+ 2 \(e\9 fi (t))\ 2 + 



(14) 



In the case of unbounded propagation in momentum space, it is the 1st term within the curly 
brackets which yields the dominant contribution at large t (the 2-nd is bounded by 1, the 3-d 
is order of square root of the 1-st.). Note that the 1-st term is the average over (3 of the (3— 
rotors kinetic energy. 
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2.3 Motion in the falling frame. 

In the case of nonzero gravity, quasi momentum is not conserved any more because the kicking 
potential in (|7]) is not periodic (unless rj is rational). However, a quasi-momentum -conserving 
evolution is restored by the substitution : 

Since 77 is the momentum gained over one period due to gravity, this amounts to measuring 
momentum in the free-falling frame (recall that the positive x— direction is that of gravitational 
acceleration). 

From (Q) it follows that 

V>/(p, t + 1) =J2 Kn z~ i¥p ~ n+tri+l)2 ^f{p -n,t) . (15) 

n 

that may be rewritten as : 

\^ f (t + 1)) = Uf(t)\l/>f(t)) , Uf(t) = e -ikco S (X) e -iZ ( P +trl +^ _ (lg) 

The operator Uf(t) describes up to a constant phase the evolution from (continuous) time 
t = (tr) + to time (t + r) + under the time-dependent Hamiltonian 



1 +00 
Hf(t) = -(P + ^t) 2 + fccos(X) «(t-«r), (17) 



t=— 00 



that describes the motion in the falling frame. This Hamiltonian is related to that in eqn.(fj) 
by the gauge transformation ^ Xx l T , The classical dynamics corresponding to ( IB ) is given by 
the area-preserving, time-dependent map in the (x,p) plane: 

Pt+i =Pt + k sin(x m ) , x t+1 = x t + r(p t +tr) + -) . (18) 
Since h = 1 in our units, the classical limit is approached as 

k — > 00 , r^O, rj — > 00 , kr — > const. 7^ , r/r — > const. 7^ . (19) 

The evolution (|TJ|) only mixes momenta which differ by integers: hence, quasi-momentum 
is conserved, so the Bloch-Wannier fibration can be fully implemented, as described in the 
previous Section. Eqn.(|T2"D now becomes: 

t-i 

!*/»(*)> = II «/»(OI*/J> > = £6/»(r) (20) 

r=0 

(all operator products ordered from right to left), where 

(n\t\m) = K n _ m , (n\K p (r)\m) = 5 nm e -^(^+n+v/2? (21) 
Consequently, in the falling frame, 

Uf,(t) = e~ ikcos( - §) e -m*+P+rit+r,/2)\ (22) 
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3 Features of Quantum Dynamics in the falling frame. 



The theoretical importance of the Kicked Rotor model lies with its asymptotic properties in 
time. Best known among these are dynamical localization, and quantum resonances. Though 
not crucial to current experiments, the long-time asymptotics is an important theoretical ques- 
tion for the present models, too. Apart from resonances, we do not attempt at a thorough 
theoretical analysis of this issue, which is likely to critically depend on the arythmetic type 
of r, rj, (3. Though some results in this respect will be presented later, we mostly focus on 
intermediate-time features of the quantum dynamics, directly connected to experimental find- 
ings. 

In this section the falling frame is constantly used, without any further reference to the 
laboratory frame. The suffix f will be omitted in order to simplify notations. 

3.1 Resonances 

3.1.1 Zero Gravity 

If the kicking period r is commensurate to Atc, that is r = Anp/q with p,q relatively prime 
integers, and in addition (3 = m/(2p) with m < 2p an integer, then the (3— Kicked Rotor 
exhibits the phenomenon of Quantum Resonance In that case, indeed, the rotor dynamics 
commutes with translations in momentum by multiples of q. This typically results in band 
(absolutely continuous) quasi-energy spectrum and ballistic spread of the rotor wave packet 
in the momentum representation. For special values of q the bands may however be flat. 
This is the case when q = 2; the ballistic spread is then only observed at (3 = 1/2, while 
at (3 = the dynamics is sharply localized in momentum ("anti- resonance"), as can be seen 
from (|T3"D, making use of the fact that e~ tnln = e~ lnln holds for integer / and n. The width of 
the quasi-energy bands rapidly decreases as the order q of the resonance increases 0. Ballistic 
motion is then observed only after quite long times. This places high-order resonances beyond 
experimental observation. Our discussion is mostly focused on the main resonances (q = 1,2). 

Since the quadratic growth of energy at resonant values of r requires a sharp value of 
quasi-momentum (hence an extended particle state in position), it cannot be observed in the 
particle wave packet long-time dynamics, not even at resonant values of r. Nevertheless, with a 
generic choice of the initial wavepacket, a resonant effect is still manifest, in the form of linear 
asymptotic growth of the particle energy. This is qualitatively understood as follows. At time 
t, rotors whose quasi-momenta lie within ~ 1/t of the resonant value(s) are still mimicking 
the resonant growth of energy oc t 2 . Assuming a smooth initial distribution of quasi-momenta, 
such rotors enter the average over j3 (|14D with a weight oc 1/t. This directly leads to the 
linear growth of E(t). The latter can also be rigorously proven, and the proportionality factor 
computed. This is done in the Appendix, under the hypothesis that the initial wave function 
of the particle is such that (ip\X 2a \ifj) < oo for some a > 1. 

3.1.2 Nonzero gravity 

The j3— rotor dynamics at resonant values of r in the presence of gravity is a nontrivial math- 
ematical problem, which may give rise to different types of quantum transport depending on 
the arythmetic type of n. 
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In the Appendix we prove that for r = 2lir ,(l integer), and for all irrational r\ in a set of full 
Lebesgue measure, the energy £(t) (in the falling frame) grows like k 2 t/4, under the hypothesis 
that the initial wave function satisfies (ip\X 2a \il)) < oo for some a > 1. 



3.2 Accelerator Modes 

3.2.1 Quantum Rotor Dynamics near Resonance. 

We shall now analyze the quantum dynamics at values of r close to the resonant values 27r/ 
(/ > integer) and for large kicking strength k. We hence assume r = lul + e with I integer 



and |e| << 1, and rescale k = k/\e\. Replacing in fl2~l|) and noting e mln = e mln we get: 

whence it follows that (apart from an irrelevant phase factor) 

U p (t) = e -^' kcos§ e -^ (i ' l) (23) 

wheref] 

I=\e\ft=-i\e\± , ?i: /3 (/,t) = ^sign( e )/ 2 + /(7r/ + r(/? + tr ? + |)), (24) 



If |e| is assigned the role of Planck's constant, then (P3j) , (P4|) is the formal quantization of either 
of the following classical (time-dependent) maps : 

It+i = h + k sm(6 t+1 ) , 9 t+1 =9 t ±I t + 7rl + t(J3 + t V + »j/2) mod(27r) (25) 

where ± has to be chosen according to the sign of e.[] The small |e| asymptotics of the quantum 
(3— rotor is thus equivalent to a quasi-classical approximation based on the "classical" dynamics 
(|25|). We emphasize that "classical" here is not related to the h — > limit but to the limit 
e — ► instead. The two limits are actually incompatible with each other except possibly when 
/ = (see (|i9|)). For the sake of clarity the term "e— classical" will be used in the following. 



3.2.2 e— Classical rotor dynamics. 

Changing variable to J t — It ± 7rZ ± r(/3 + rjt + rj/2) removes the explicit time dependence of 
the maps fl25D , yielding: 

Jt+i = J t + ksm(9 t+1 )±rr] , 9 t+1 = 9 t ±J t . (26) 

These area-preserving maps are 27r- periodic in J and 6, so they define smooth, area-preserving 
maps of the 2-torus parametrized by J = J mod(27r), d = 6 mod(27r). Such "toral maps" are 
conjugated to each other under J — > 2n — J mod(27r), •& — > $ + 7r mod(27r). They differ from 
the Standard Map by the constant drift tjt. 

4 The integer time variable t is fixed on both sides, and the 2nd exponential on the right hand side is that of 
a constant operator. 

5 The double sign would not appear if e were used as the "Planck constant", rather than |e|. In this way the 
sign of I would be reversed with respect to that of physical momentum whenever e < 0. 
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Figure 1: Phase portraits for the map ( f26| ) on the 2-torus, with r = 5.86, rj/r = 0.01579, e = —0.423. 
Here the torus is mapped onto [—tt, tt) x [— tt, tt). One stable fixed point with j = exists for 0.542 < 
k < 4.037. Left: k = 1.329, the stable fixed point is at J* = 0, 9* = 0.42. Center: k = 4.232, a stable 
period-2, j = orbit is visible. Right: k = 4.494, the period-2 orbit has left room to a stable period-6 
orbit. The lower left part is magnified in the inset. 



Let j7~o, $o be a period-p fixed point of either toral map. Iterating (|^) with J = Jq , $o = 
i9 , at t = p one obtains: 

J p = J + 2vrj , e p = 6 + 2vm (27) 

for some integers j , n. The points {Jt^t) with < t < p — 1 are period-p fixed points 
themselves, and define a period-p periodic orbit of the toral map, which is primitive if all such 
points are distinct. 

Period-1 fixed points are given on the 2-torus by j7o = 0, i?o = 9j, where 

M e,) = (28) 

k 

and j is any integer such that: 

-k ± tt] < 2ix) < k ± tt] . (29) 



No such point exists if k < rrj < tt; two (at least) exist (that is, (|28|) is solvable for at least one 
value of the integer j) whenever k > it. Two period-1 fixed points with j = exist whenever 
k > tt]. In order for the period-1 fixed points ( p8|) to be stable it is required that: 

-4 < ±fccos($ ) < (30) 
From (E3), (|30D it follows that for any integer j each map (BB) has exactly one stable period-1 



fixed point on the 2-torus, given by (|28|) if, and only if, 

< k < , k$£ = \2tt) T r V \ , kSS = v / 16 + (27rj T rr ? )2 . (31) 
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Figure 2: Same as Fig.[j], for rj/r = 0.01579, k = 0.87r(47r — r), and different values of r. Left: at 
r = 10.996, one stable fixed point exists with period 1, j = 0. On decreasing r it turns unstable at 
r rj 10.813, generating a period-2 stable orbit. At r « 10.801, one stable orbit with period 1, j = — 1 
appears. Both orbits are shown in the phase portrait on the right, drawn at r = 10.744. 



At k = k\^ ax such fixed points turn unstable and bifurcations occur. This is shown in Fig.|T] 
for a case with j = 0. At k = k^~ x = 4.037 a stable period-2 orbit appears. This becomes in 
turn unstable at ~ 4.490, and a stable period-6 orbit is left . The size of the islands rapidly 
decreases through the sequence of bifurcations. At k = 5 no significant stability islands are 
any more visible; at k = 5.741 and at k = 6.825 stable period-1 points with j = —1 and j = 1 
respectively appear. The rise and fall of these, and of subsequent higher-j period-1 points as 
well, are ruled by ([31]). Further examples of period-1 fixed points, and bifurcations thereof, are 
illustrated in Fig]^. 

Examples of primary (that is, not born of period-1 points by bifurcation) higher-period sta- 
ble orbits are shown in Fig.§. Generally speaking, the presence of two independent parameters 
(k and rrf) in the maps produces a remarkable variety of stable periodic orbits of higher periods, 
depending on the parameter values in complicated ways. Figs.0, [| and |3| provide but a partial 
view of such complexity. They were singled out because they pertain to experimentally rele- 
vant parameter ranges; see the discussion in sect. (3. 2. 6). In particular, the value rj/r — 0.01579 
constantly used in this paper is that of experiments in |T7|, [18], [L9| . 



3.2.3 e— Classical Accelerator Modes. 

Primitive periodic orbits of the toral maps yield families of accelerator orbits of the original 
dynamics (25), marked by linear average growth of momentum with time: 

9pt = 8o = #o mod(27r) , I tp = 7 + apt , (32) 
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Figure 3: Phase portraits for the e-classical dynamics on the torus, showing a (5, —2) periodic orbit 
at r = 12.472, k = 0.8vr(47r — t), tj = 0.01579r (left) and a (10,1) periodic orbit at r = 6.31, 
k = 0.8vr(r - 2vr), n = 0.01579r (right). 



where 




(33) 



with m any integer, and Jq^q a period-p fixed point. If such accelerator orbits are stable, 
then they are surrounded by islands of positive measure in phase space, also leading to ballistic 
(linear) average growth of momentum in time. These are named accelerator modes. 

We shall classify accelerator modes according to their order p and jumping index ). By a 
(p, j)-accelerator mode we shall mean a mode, whose order and jumping index are given by the 
integers p,j respectively. 

3.2.4 Quantum Accelerator Modes in Rotor Dynamics. 

Initial physical momenta n = |e| _1 / for e— classical accelerator modes are obtained from ( |3"B| ) 
for any < (3 < 1 : 



If the stable islands associated with e— classical accelerator modes have a large area compared 
to |e|, then they support a large number of quantum states. Thus they may trap some of the 
rotor's wave packet and give rise to quantum accelerator modes traveling in physical momentum 
space with speed ~ a/|e| = —rrj/e + 27rj/(p|e|). In order that such modes may be observed, 
the phase space distribution associated with the initial rotor state must significantly overlap 
the islands. Even in that case the modes will eventually decay due to quantum tunneling out 
of the classical islands. 




(34) 
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Figure 4: Quantum phase space evolution of the (3— rotor with /? = 0.218769, k = 1.329, r = 5.86, 
77 = 0.093 initially prepared in the coherent state centered at the e— classical (1,0) accelerator mode 
no = —3.754, 9q = 0.420. Contour plots of the Husimi functions of the rotor are shown at times 
t =0,2,4,8,16. The black spots in the centers of the contours are an ensemble of classical phase points, 
initially distributed in a square of area h = 1 centered at the mode. They evolve according to the 
e— classical dynamics (^5|). 



This picture is confirmed by numerical simulations. Fig.f| shows the quantum phase-space 
evolution of a (3— rotor started in a coherent state centered at the position of the (1, 0)- acceler- 
ator mode generated by the fixed point shown in Fig.[l] (left). The Husimi functions computed 
at subsequent times closely follow the motion of the e— classical mode. 

In Fig.[5| we compare quantum and e-classical results after 30 kicks, and different values of 
r at fixed k and rj/r. In some of the examined r range, two different e-classical modes are 
simultaneously present, with different signs of the acceleration. In order to single out one of 
them, we have plotted the quantities E±, which, in the classical case, are equal to the average 
energy after 30 kicks, computed over those rotors in the ensemble which, at the given time, 
have a positive (respectively, negative) momentum. In the quantum case, they are equal to the 
energy, computed in the state obtained by projecting the rotor state over the positive (respe., 
negative) momenta. The main peak in the left-hand plot is due to a (1, 0) mode whose interval 
of existence and stability is, according to (|3"TD, 1.043 < t/(2tt) < 1.261. Both the classical 
and the quantum data sharply arise at the onset of the mode. The rise of the quantum data 
is actually milder, because the e-classical island has to grow beyond a size ~ e in order to 
be quantally significant (note that |e| increases on increasing r). As the stability border is 
approached, the island starts shrinking again, the classical and the quantum mode become less 
and less effective; the latter somewhat faster, for the just mentioned reasons. At the stability 
border, a (2,0) classical mode is originated by bifurcation, which is in fact signalled by a tiny 
peak in the classical data; however, no similar quantum structure is observed. The small peak 
on the leftmost part is an accelerator mode itself, see Sect. (3.2.6). 
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Figure 5: Full dots: Quadratic spread E±(i) of the quantum /3-rotor over positive (right) and negative 
(left) momenta after 30 kicks, vs r, with t]/t = 0.01579, k = 0.87T, f3 = lir/r — rj/2, I = 1 (left) and 
I = 2 (right). The rotor was started at n = 0. Full lines: same for an ensemble of 5 x 10 6 classical 
rotors evolving according to (25), started with n = and 9 uniformly distributed in [0, 2ir]. The exact 
meaning of E±(t) is explained in the text. At values of t/2tt ~ 1.72 the classical phase portraits are 
as shown in Fig. 2. 



At not really small values of the "Planck's constant" e, strong e-classical modes still leave 
quantal signatures. The quantum-classical correspondence turns however loose, in interesting 
ways. This is illustrated in the right hand part of FigJ^. Now |e| decreases from left to right. 
The rightmost peak corresponds to a (3, —1) mode (see Sect. 3.2.6), faithfully reproduced by 
quantum data. The main classical peak is the (1, 0) mode, active for 1.721 < t/(2it) < 1.862. A 
significant quantum mode is also observed, with remarkable differences, however. In particular, 
the quantum mode has its maximum at the very point where the (1, 0) e-classical mode turns 
unstable, giving birth to a (2,0) mode, as demonstrated in Fig. 2 (right). Our explanation of 
this curious behaviour is as follows. When the stability island near a stable fixed points breaks 
at a bifurcation point, some remnants of its KAM structure nevertheless survive, in the form 
of broken tori (cantori). These provide but a partial barrier to classical motion, and allow 
for nonzero phase-area flux. If this flux is small compared to e, then the cantorus quantally 
acts as if it were unbroken [22] . Then it may give rise to a quantum accelerator mode, much 



more effectively than one might guess looking at the small area of stability of the bifurcated 
higher-period orbits. 

Though e— classical modes exist at any value of f5, their location in the /3-rotor's phase-space 
changes with f3 . Hence their impact on the quantum evolution of a rotor state is enhanced at 
those values of (3 which afford maximal overlap of the mode with the given state. In particular, 
if t]t < 2tt, and the (3— rotor is initially set in the n = momentum eigenstate, then the 
quantum (1, 0)-accelerator modes are especially pronounced when (3 = ln/r — rj/2. Note that 
(3 was not set to this value in the case of of Fig.f|. 

In the e— semiclassical regime, accelerator modes exponentially decay in time due to quan- 
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Figure 6: Exponential decay of accelerator modes. Left: semilogarithmic plot of the probability P$ n (t) 
in a moving momentum window of width 8n = 15 centered at the e— classical (1, 0) mode , versus time 
t, for k = 2.751, r = 5.8, e = -0.483, rj = 0.092, (3 = tx/t - r?/2. The asymptotic decay of P 5n (t) 
is exponential with decay rate j € ~ 0.32 x 10 -4 . Right: semilogarithmic plot of 7 e vs l/(|e|). The 
e— classical variables 51 = \e\5n, k = \e\k were held fixed at 7.540, 1.329 respectively while changing e. 
The fitting line corresponds to the law 7 e ~ 135 x exp(— 7.35/|e|). 

turn tunneling out of the classical islands, with decay rate 7 e oc exp —(const. /|e|). The decay 
of quantum accelerator modes is numerically demonstrated in Fig.|6]. 



3.2.5 Quantum accelerator modes in Particle Dynamics. 

Quantum Accelerator modes arise in the particle dynamics as well, just because such dynamics 
comes of a quantal superposition of j3— rotors. We shall presently discuss the small-e asymptotics 
of the particle dynamics. This will at once provide an alternative derivation of the e-classical 
rotor dynamics, and a means of translating into particle dynamics the results established in 
the previous sections. 

Let us consider the propagator from state \p) to state \x') from (discrete) time t to time 
t+1, for the particle dynamics (]TB|). Let p = n + f3 as usual, and x' = 2m' it + 6' with m' integer 
and < 9' < 2tt. Then 



(x'\U(t)\p) 



2,71 



exp(i0(/?,0',t)) 



exp ( — ik cos(6 |/ ) + inO' + 2im'n/3 en 2 — innl — inr((3 +tij-\ — ) ) ,(35) 



where e 



i-nln 



was used, and 



(36) 
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Next we introduce e-classical scaled variables I = n\e\, k = k\e\, V = —2irm'\e\, to be kept 



constant in the e-classical limit. Then (|35|) rewrites as: 



2tt 



where: 

T{9', L', ft, I, t) = 19' -L'ft-k cos(0') - ^sign(e)/ 2 - I {trl + r(ft + tr) + r]/2)) . (38) 

Considering |e| as the Planck's constant, and J, L' as canonical momenta respectively conjugated 
to 9, ft', the function T is a generating function for the canonical transformation (9, 1, ft, L) — > 
(9', V, ft', If) given by: 

' - 

r = w = I + ' ksi " lt>,) 

dT 

9 = — = 9'- sign(e) I - irl - r(ft + trj + rj/2) . (39) 

The second exponential in ( |3"7D is thus (apart from a constant prefactor) the e-semiclassical 
propagator associated with the e-classical map (|39|). The 3d and the 4th equation are just the 
e-classical /3-rotor dynamics. The first equation says that quasi-momentum is conserved, and 
the second yields the complete revolutions accumulated by the /3-rotor from time t to t + 1; 
this quantity is formally conjugated to quasi-momentum. 

However, (|37D has the additional phase factor e 1 ^ , which is not e-semiclassical, because is 
not scaled by |e| _1 . While this factor is irrelevant for the fixed ft dynamics (and was in fact 
disregarded in sect. (3.2.1)), it cannot be neglected when studying the particle wavepacket 
dynamics, which requires integration over all values of ft. Such integration causes e-classical 
trajectories of different /5-rotors to interfere. This interference is ruled by the true Planck's 
constant fi = 1, and cannot be suppressed by the e — > limit. Unlike the /3-rotor dynamics, 
the particle dynamics does not become "classical" in the e — > limit. 

As long as ft is fixed, the maps ([39]) yield, in the physical variables p, x and at \p\, \x\ » 1: 

p t = \e\- 1 I t + ft ~ Po + TT7] \ 



VP / M 

,-1 r f 2 ^ \ Tt(t-1) . s 

x t ~ -|e| 1 L t ~ Xo + n Tt+ ^— Trr]J K — (40) 

for a (p,j)- accelerator mode (|32|) started at Xq, po. 

Fig.0 shows the quantum evolution of a particle started at t = in the coherent state 
centered at the (1, 0)-accelerator mode p = 0, x = 0.42. The phase-space distribution splits, 
the righmost part of it moving with constant acceleration —rrj/e. This is the effect of the 
accelerator mode. 
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Figure 7: Quantum phase space evolution for a particle initially prepared in the coherent state 
centered at po = 0, xq = 0.42, with k = tt, t = 5.86, r\ = 0.093. The computation implements the 
Bloch-Wannier fibration discretized over 512 values of quasi-momentum. Shades-of-grey plots of the 
Husimi function of the particle at times t = 2,8, 16 are shown. 

A conceptually simpler situation is met when the initial state of the particle is an incoherent 
mixture of plane waves, for in that case no interference occurs between different j3— rotors. Let 
the initial particle state be described in the falling frame by the statistical operator: 



P(0) 



dp f(p)\p)(p\ , f(p) > , dp f(p) = 1 . 



(41) 



Each plane wave has a well-defined quasi-momentum, so it is equivalent to a unique (3— rotor 
in the angular momentum eigenstate specified by the integer part of the momentum of the 
wave. Therefore, the statistical ensemble (f41~[ ) is equivalent to a statistical ensemble of (3- 
rotors. At any given quasi-momentum {3, the state of the rotor is described by the statistical 
operator p^ = (P(/3)) _1 J2 n f( n + P)\ n )( n \i the distribution of quasi-momenta is further given 
by P{(3) = ^2 n f(n + (3). The momentum distribution for the particle is given at time t by 



f(p,t)=P((3)(n\pp(t)\n) 



(42) 



where (3 = {p}, n = [p], and f>p(t) evolves according to the /3-rotor dynamics ( p0|) . The 
distribution in momentum is then e-quasi classically that of an ensemble of classical rotors 
evolving according to (|25|). 



3.2.6 Spectroscopy of accelerator modes. 



In the experiments described in refs. |I7|, |18| , |19[, the initial state of the falling atoms is sat- 
isfactorily described, according to the same references, by (|4"ID , with f(p) a Gaussian of rms 
deviation ~ 2.55 (in our units) centered at p — 0. In Figs. |^ and |9] we show numerical results 
obtained with the same choice of the initial state. Such results provide further evidence of 
accelerator modes, including higher order ones, which were not previously identified. The Fig- 
ures were produced by computing the evolution of an ensemble of 50 rotors with the mentioned 
initial distribution over 60 kicks, for different values of the period r near the main resonances 
r = 2n, r = 47r, and k = 0.87T. As rj/r was kept fixed at the physical value 0.01579, r\ also 
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Figure 8: Momentum distribution in the falling frame after 60 kicks for different values of the kicking 
period r near r = 2tt, and for k = 0.87T, r/ = 0.01579r. Note the negative sign of p. Darker regions 
correspond to higher probability. The initial state is a mixture of 50 plane waves sampled from a 
gaussian distribution of momenta. Full lines are the theoretical curves (|43|); they are drawn piecewise, 
to avoid hiding the actual structures to which they correspond. Their order and jumping index are 
indicated by the arrows. The inset shows data at t = 400. 
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Figure 9: Same as Fig.@, for r near 4-7T. The inset shows data at time t = 400. 



varied with r. Fig|| is analogous to the experimentally obtained Fig.2 in ref. [|17| (with some 
differences in units and in parameter values, though). 

The final distribution of momenta is represented by shades of grey in the (p, r) plane, darker 
zones corresponding to higher probability. The hyperbolic-like structures near r = 2tt,Att are 
signatures of quantum accelerator modes. At any r where they are visible, they are in fact 
located at the momentum values reached at t = 60 by certain accelerator modes, started at t — 
with n ~ 0. Such values are theoretically predicted by eqs. fl32|) . The e-classical (pj)-mode 
started at t — with Iq = no lei is located at time t at the momentum: 



n 



Til 

n ~ t— + t 



2tt] 



(43) 



Replacing e = r — 2nl (I — 1, 2 for Fig.|| and Fig.|9] respectively), and r\ = 0.01579r, one obtains 
a curve n = n{r) for any chosen time t and for any chosen no, p, j. Such curve is then observable, 
if a mode with the chosen p,j exists, which significantly overlaps the initial distribution. The 
narrow distribution of initial momenta legitimates the choice no = 0; although other small 
values of no are involved, they just result in a thickening of the curves. 

Well-marked (1,0) modes are observed in both Figures. The intervals of existence and 
stability of the e-classical (1,0) modes predicted by ( |31~D are 0.745 < t/(2tc) < 0.963 and 
1.043 < t/(2tt) < 1.261 for the case of Fig.§. The (1,0) mode in Fig.9 is the same as in Fig.5; 
in most of its range the e-classical structure is like the one shown in Figs|| and [5[ It has both a 
(1, —1) stable orbit and a (2,0) stable orbit originated by bifurcation. The (1, —1) theoretical 
curve mostly lies at negative values of p off the scale of the figure, and no significant trace of it 
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was detected in our quantum computation, indicating that the island was too tiny compared 
to the relatively large values of |e|. We therefore explain the large (1, 0) structure obaserved by 
the same quantal mechanism discussed in sect. (3.2.4), namely localization by cantori. 

Higher-order modes (10, ±1) and (5, ±2) are observed close to r = 2tt and r = An respec- 
tively. The corresponding e-classical modes are shown in Fig|| The small, yet well marked, 
structure produced by the (10, 1) mode near 2n is also visible in experimental data in ||17|| . The 
correspondence of the (5, 2)- and (10, l)-modes with the curves fl43|) is remarkably evident at 
longer times, see the insets in Figs|| and |9|. 

In Figs.[8] and |9| other higher-order modes are visible, too. These were identified by first 
fitting the observed structure with a curve and then checking that the e-classical phase 
space really displays, in the relevant r range, a stable orbit with the found p, j. A few of these 
left but a dim trace in Figs.[5] and |^. This may be due on one hand to mutual interference of 
different modes when they coexist in the same r-range, and on the other hand to the small 
number of rotors used in the computation. For such reasons, Figs|| and | do not probably 
account for all the accelerator modes which are excited in their respective parameter ranges. 

Generally speaking, the family of higher-order modes that are potentially observable in 
Figures like [| and |9] (and in experiments as well) is probably much richer than shown here. 
These might be exposed by varying parameters, and also by a fine scan of smaller e ranges. 
It looks likely that momentum distributions at relatively short times, of the type shown in 
Figs.|8| and can be altogether described by accelerator mode spectroscopy, i.e identification 
of accelerator modes and analysis of their mutual interference. Such a systematic analysis is 
beyond the scope of this work. 



4 Concluding remarks. 

1) The above theory of accelerator modes hinges on reduction to independent kicked rotors 
models, wherein such modes admit a natural interpretation in terms of classical trajectories. 
It is the conservation of quasi-momentum that allows for such reduction. Accelerator modes 
in particle dynamics are a purely quantal effect, in fact a remarkable manifestation of the 
conservation of quasi-momentum (in the falling frame). 

2) The fact that the intermediate-time dynamics is dominated by a discrete set of modes, which 
exponentially decay in time, bears a distinct resemblance to the Wannier-Stark problem of a 
Bloch particle in a constant field [23]. How far this analogy carries is, in our opinion, a very 



interesting question. 

3) The time- dependent variants of the Kicked Rotor model examined in this paper also raise 
other important theoretical questions, which were not addressed in this paper. These are 
about the long-time asymptotics of the dynamics, and the localization-delocalization issue in 
particular. 

4) In the experiments [17], [18, [19[] the initially prepared states form an incoherent superposition 



of momentum eigenstates, and of quasimomentum states. Therefore the results of the of the 
calculations for the /3-rotors explain the experiments. The particle nature (compared to the 
rotor) is important only for initial coherent superpositions of momentum eigenstates. It will of 
great interest if this point is explored experimentally. 
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5 Appendix: Resonances in the presence of Gravity 

In this appendix we assume r = 2irl, with I ^ integer. We denote ip(x) = (x\ip) the wave 
function of the particle in the position representation, and S(t) the kinetic energy at time t (in 
the falling frame). 

Proposition: Let J dx \x\ 2a \il){x)\ 2 < +00 for some a > 1. Then: 
(I) If ri = Q : 

k 2 Dt 1 
S(t) = £(0) + — + 0(t X ) , X = max{2 - a, -} 



where: 



n=0 ^° 



I 2 , Pn = j + \ mod(l) . (44) 



(II) If rj ^ satisfies a Diophantine condition, that is, there are constants c, 7 > so that for 
all integer p, q : 



p 

T}-~ 

Q 



> cq- 2 ^ (45) 



then: 



k 2 t I 
£{t) = £(0) + — + 0(t a ) , a = max{2 - a + 7, -) (46) 



Remarks: 

1. The quasi- momenta /3 n 's in part (I) are exactly those yielding quadratic growth of the 
j3— rotor energy. 

2. Part (II) ensures that, for all rj in a set of full measure, the energy grows diffusively with 
coefficient k 2 /4. 

Proof: Numerical constants will share the common notation C whenever their exact value is 



irrelevant. From (2~H) 



The explicit form of the phase (J) is not important for our present purposes. Further, 

(e\y p (t+i)) =(e\fcK f) (t)\Vf i (t)) = eWi,i) e -*™(«)(j_fl_w|f^)) , 

where: 

a = vr/(l + 2(3 + rj) , b = 2ixlr] . (47) 

It follows that: 

(6\v p {t)) = e -imi,t) e -ikF(8,p,t)( e _ at _ K t _ ^1^(0^ 
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where: 



t-i 



F(6,P,t) = ^cos{9 -ra-rbt + r 2 b/2 + rb/2) 



r=0 



We now use eqn.(|14]). As already remarked, the dominant contribution is given by the 1st term 
on the rhs, corrections being on the order of square root of that term. We hence restrict to 
that term. After substituting the above equations, it takes the form: 



dp 



2k 



de 



£ {t)+S 1 {t)+£ 2 {t), 



(48) 



having denoted 



1 r 1 

2 



^dpj^de\g\9M 



£o(t) = 

= h i!ldpjl*d9F">{e,p,t)\g{9,t)\\ 

£ 2 {t) = jfeK/o 1 dp J 2 " d9 iF'(9,p,t)g*(9,t)g'(9,t) 
g(9,t) = (9 -ta-bt 2 /2 + bt/2\^ f3 (0)) 

(primes denote differentiation with respect to 9). An obvious shift in 9 shows that 

£o(t) = const = £q(0). 

Moreover, from the Cauchy-Schwarz inequality it follows that: 



(49) 



|£ 2 (t)| < k(£ 



27T \ 1/2 

dp I d9F' 2 (9,P,t)\g(9,t)\ 2 



1 p2ir 

dp / d9\g'{9,t)f 
Jo 



1/2 



2E 1 (t) l/2 £M 



/2 



(50) 



The dominant contribution to £{t) is thus given by £\(t), so £ , £ 2 will not be considered until 
the end of the proof. After an obvious change of variables, 

2 



2 



l di3 l <Hs 



sm(9 + ra + br 2 /2-br/2) | (9\^ p (0)) | s 



(51) 



We rewrite the square of the sum as a double sum, then apply standard trigonometric formulae, 
and finally replace a, b by (^7|) , leading to : 



k 2 f 1 f 2n 



with 



^(0))| 2 (^(/^-^OW)) 



1 £ (_)(r+-)i cos (2vr//?(r - s) + 7ilr](r 2 - s 2 )) , 



(52) 



r,s=l 



Ai(P,t) = 

A 2 (p, 9,t) = |E (-) {r+s)l cos (29 + 2vr/3(r + s) + nlr](r 2 + s 2 )) 

r,s=l 

Now we expand |(^|^ r j g(0))| 2 in Fourier series: 

| (01*0(0)) | 2 = c(M,N) e 2niMp e l 



AN 6 



(53) 



M,N 
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Replacing in (|52|) we obtain 

1 



E x (t) = -nk^B, - B 2 ) , (54) 



where: 



B 1 = c{l{r-s),0){-Y r - s)l e~ i7Tlr ^ 2 - s2) (55) 

r,s=l 

t-l mm(t,j+t) 
= c{l],0){-) lj e lnlvj2 e" 2 ™^' , 

j=l-t r=max(i,j+i) 
t 

B 2 = c ( / ( r + s )' 2)(-)^ +s ) e -^ (r2+s2) 

r,s=l 

2* min(ij-i) 
= J2c{l],2){-) lj e- inlrij2 e - 2 ™ lrir{r - j) . (56) 

i=2 r=max(ij-t) 

With the help of the Lemma proven in the end of this section, B 2 is bounded by : 

1+t 2t 2t 2t 

J2\ c (ij,2)\(j-i)+ \c(ij,2)\(2t- 3 + i)<cJ2(j-i)r a + ctY,r a = o(?- a ). 

j=2 j=t+2 j=2 j=t+2 

Here and in the following 0(t x ) has to be read as O(logt), 0(1) whenever x = 0, x < 
respectively. In order to estimate Bi we distinguish two cases; 

Case I: 1] = 0. 

t—l oo 

B 1 = J2c(lj,0)(-) lj (t-\j\) = tY,(-) lJ c(lj,0) + O(t^) 

j=l-t j=~oo 

= ^ + °(^ Q ) (57) 



where 



D = 2n J2 



j=-oo 
l—l oo 



n=0 j'=—oo 



-J2 d6\(9\^^ , (3 n = j + - mod (1). (Ot 

n=0 ^° 



We next substitute fl57|) and (^) in (|54|) , and then in (f48|). Recalling (f|^),(^(]), and the remark 
preceding (f48l) leads to 
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Case II: rj a Diophantine irrational. The sum (|55D is written in the form 

B x = c(0,0)t + S 
where S is the contribution of all j =^ terms. It can be bounded as 

~ 2 ^' C ^ 7 '°' ) 'l sin ( 7r %)l 

t-i 

< CY,mM3 l+1 = 0(?- a ^). (59) 

where (|45|) was used. The normalization of the wavefunction implies c(0, 0) = l/(27r). Substi- 
tuting in (|>|) leads to fllB]) after the same concluding steps as in Case I above. This completes 
the proof. □ 

Lemma. Under the hypotheses of the Proposition, and with c(M, N) defined as in (\5^), 
\c(M,N)\ = 0(\M\- a ) as M -► oo. 

Proof. From the Bloch-Wannier fibration (|) it follows that, if \M\ > 1, then: 

\c(M,N)\ = — I dxi)*{x + M-K)i)(x- M-k) e~ iNx 
2n J 

< C\M\~ a J dx\i)*{x + M-K)i){x-M-K)\Q{x,M) (60) 

where Q(x, M) = (1 + (x + Mvr) 2 ) Q / 2 (l + (x - Mir) 2 ) a / 2 > (An 2 M 2 ) a / 2 . The Cauchy-Schwarz 
inequality then yields : 

\c{M,N)\ < C\M\~ a J dx \ij(x)\ 2 (1 + x 2 ) a < oo , 

as the convergence of the integral was assumed in the Proposition. □ 
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